Quantal molecular description and universal aspects of the spectra of bosons and 

fermions in the lowest Landau level 



Constantine Yannouleafl and Uzi Landmar0 

School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332-0430 

(Dated: 29 November 2009) 

Through the introduction of a class of trial wave functions portraying combined rotations and 
vibrations of molecules formed through particle localization in concentric polygonal rings, a cor- 
related basis is constructed that spans the translationally invariant part of the lower-Landau-level 
(LLL) spectra. These trial functions, referred to as rovibrational molecular (RVM) functions, gen- 
eralize our previous work which focused exclusively on electronic cusp states, describing them as 
pure vibrationless rotations. From a computational viewpoint, the RVM correlated basis enables 
controlled and systematic improvements of the original strongly-correlated variational wave func- 
tion. Conceptually, it provides the basis for the development of a quantal molecular description for 
the full LLL spectra. This quantal molecular description is universal, being valid for both bosons 
and fermions, for both the yrast and excited states of the LLL spectra, and for both low and high 
angular momenta. Furthermore, it follows that all other translationally invariant trial functions 
(e.g., the Jastrow-Laughlin, compact composite- fermion, or Moore- Read functions) are reducible to 
a description in terms of an excited rotating/vibrating quantal molecule. 
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I. INTRODUCTION 



Motivation 



Following the discovery [l[ of the fractional quantum 
Hall effect (FQHE) in two-dimensional (2D) semiconduc- 
tor heterostructures under high magnetic fields (B) in 
the 1980's, the description of strongly correlated elec- 
trons in the lowest Landau level (LLL) developed into 
a major branch of theoretical condensed matter physics 
[2l4l9l|. Early on, it was realized that the essential many- 
body physics in the LLL could be captured through trial 
wave functions. Prominent examples are the Jastrow- 
type Laughlin (JL) [3], composite fermion (CF) Q, and 
Moore and Read's (MR) J] Pfafhan functions, represent- 
ing quantum- liquid states 0- In the last ten years, the 
field of semiconductor quantum dots [l^ helped to focus 
attention on finite systems with a small number (A'') of 
electrons. Theoretical investigations of such finite sys- 
tems led to the introduction of "crystalline" -type LLL 
trial functions referred to as rotating electron molecules 
[12, (REMs). In particular in their intrinsic frame 
of reference, the REMs describe electrons localized at 
the apexes of concentric polygonal-ring configurations 
(ni, n2, n^), where X]q=i = ^ ^^'^ ^hc num- 
ber of concentric rings. 

More recently, the emerging field of graphene quantum 
dots [13, Hlj , and the burgeoning field of rapidly rotat- 
ing trapped ultracold neutral gases [23 - l3^ have gener- 
ated significant interest pertaining to strongly correlated 
states in the lowest Landau level. Furthermore, it is an- 



ticipated that small (and/or mesoscopic) assemblies of ul- 
tracold bosonic atoms will become technically available in 
the near future |32h35| | and that they will provide an ex- 



cellent vehicle 




for experimentally reaching 



'Electronic address: 
1 Electronic address: 



Constantine.YannouleeisQphysics.gatech.edu I 



Uzi . Landman@phy sics . gatech. eduf " 



exotic phases and for testing the rich variety of proposed 
LLL trial wave functions. 

Despite the rich literature and unabated theoretical 
interest, a unifying physical (as well as mathematical) 
description of the full LLL spectra (including both yrast 
[36j and all excited states), however, is still missing. In 
this paper, a universal theory for the LLL spectra of a 
finite number of particles valid for both statistics (i.e., 
for both bosons and fermions) is introduced. The LLL 
spectra are shown to be associated with fully quantal (STj 
and strongly correlated ro- vibrational molecular (RVM) 
states, i.e., with (analytic) trial functions describing vi- 
brational excitations relative to the set of the special 
yrast states known as cusp states. 

The cusp states exhibit enhanced stability and magic 
angular momenta (see below), and as such they have at- 
tracted considerable attention. However, the cusp states 
represent only a small fraction of the LLL spectrum. The 
molecular trial functions associated with them are purely 
rotational (i.e., vibrationless) and were introduced for the 
case of electrons in Ref. [121 under the name rotating 
electron molecules (REMs). The corresponding purely 
rotational bosonic analytic trial functions for cusp states 
[called rotating boson molecules (RBMs)] are introduced 
in this paper; see Section HL^ More importantly, this pa- 
per shows that the quantal molecular description can be 
extended to all other LLL states (beyond the special cusp 
states) by introducing (see Section Hi Dl) analytic expres- 
sions for trial functions representing ro-vibrational exci- 
tations of both REMs and RBMs. These ro-vibrational 
trial functions include the REM or REM expressions as 
a special case, and they will be referred to in general as 
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RVM trial functions. 

It is remarkable that the numerical results of the 
present theory were found in all tested cases to be 
amenable (if so desired) to an agreement within machine 
precision with exact-diagonalization (EXD) results, in- 
cluding energies, wave functions, and overlaps. This nu- 
merical behavior points toward a deeper mathematical 
finding, i.e., that the RVM trial functions for both statis- 
tics provide a complete and correlated basis (see below) 
that spans the translationally invariant (TI) subspace 
Q of the LLL spectrum. An uncorrelated basis, with- 
out physical meaning, built out of products of elemen- 
tary symmetric polynomials is also known to span the 
(bosonic) TI subspace [ssj . 

For the sake of clarity, we comment here on the use of 
the terms "correlated functions" and/or "correlated ba- 
sis." Indeed, the exact many-body eigenstates are cus- 
tomarily called correlated when interactions play a domi- 
nant role. Consequently, a basis is called correlated when 
its members incorporate/ anticipate effects of the strong 
two-body interaction a priori (before the explicit use of 
the two-body interaction in an exact diagonalization). In 
this respect, Jastrow-type basis wavefunctions (e.g., the 
Feenberg-Clark method of correlated-basis functions [39- 
and/or the composite-fermion basis [18, 19, 42, 43|) 
are described as correlated, since the Jastrow factors in- 
corporate the effect of a strong two-body repulsion in 
keeping the interacting particles apart on the average. 
Our RVM basis is referred to as correlated since, in addi- 
tion to keeping the interacting particles away from each 
other, the RVM functions incorporate the strong-two- 
body-repulsion effect of particle localization in concentric 
polygonal rings and formation of Wigner molecules; this 
localization effect has been repeatedly demonstrated via 
EXD calculations in the past decade (see, e.g., the review 
in Ref. [l^ and references therein). In this spirit, we de- 
scribe the basis of elementary symmetric polynomials as 
"uncorrelated," since the elementary symmetric polyno- 
mials do not incorporate/anticipate this dominant effect 
of a strong two-body repulsion, i.e., that of keeping the 
interacting particles apart. 

We are unaware of any other strongly-correlated func- 
tions which span the TI subspace. Indeed, although 
the Jastrow-Laughlin function (used for describing yrast 
states) is translationally invariant, its quasi- hole and 
quasi-electron excitations are not ^ . Similarly, the com- 
pact composite-fermion trial functions are translationally 
invariant [s^l , but the CF excitations which are needed to 
complete the CF basis are not [11 [H, IH, El] . The short- 
coming of the above well known correlated LLL theo- 
ries to satisfy fundamental symmetries of the many-body 
Hamiltonian represents an unsatisfactory state of affairs, 
and the present paper provides a remedy to this effect. 
In this context, we note that although the Moore- Read 
functions 0] are also translationally invariant, they ad- 
dress only certain specific LLL states and they do not 
form a basis spanning the TI subspace. 

Our introduction of a correlated basis that spans the 
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FIG. 1; (Color online) LLL spectra for TV = 3 scalar bosons 
calculated using exact diagonalization. Only the Hamiltonian 
term containing the two-body repulsive contact interaction, 
g&izi — Zj), [see Eq. ((3|] was considered in the exact diag- 
onalization. The gray solid dots (marked by arrows; green 
online) denote the translationally invariant states. The dark 
solid dots are the spurious states (see text). The dashed line 
denotes the yrast band, while the cusp states are marked by 
a circle. Energies in units of (//(vrA^). The number of transla- 
tionally invariant states is much smaller than the total number 
of LLL states. 



TI subspace is of importance in the following two ways: 
(1) From a practical (and calculational) viewpoint, one 
can perform controlled and systematic stepwise improve- 
ments of the original strongly-correlated variational wave 
function, e.g., the pure REM or RBM. (For detailed il- 
lustrative examples of the rapid-convergence properties 
of the RVM basis, see the Appendix.) This calculational 
viewpoint was also the motivation behind the introduc- 
tion of other correlated bases in many-body physics; see, 
e.g., the treatment of quantum liquids and nuclear mat- 
ter in Refs. [39-411 and the composite-fermion correlated 
basis in RefsHll El, El IH • (2) Conceptually, it guar- 



antees that the properties of the RVM functions, and 
in particular the molecular point-group symmetries, are 
irrevocably incorporated in the properties of the exact 
LLL wave functions. Furthermore, it follows that all 
other translationally invariant trial functions (e.g., the 
JL, compact CF, or Moore- Read functions), are reducible 
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FIG. 2: (Color online) LLL spectra for A'' = 4 spin-polarized 
electrons calculated using exact diagonalization. Only the 
Hamiltonian term containing the two-body Coulomb interac- 
tion, e^/(K|zi — Zj\), [see Eq. ([3])] was considered in the ex- 
act diagonalization. The gray solid dots (marked by arrows; 
green online) denote the translationally invariant states. The 
dark solid dots are the spurious states (see text). The dashed 
line denotes the yrast band, while the cusp states are marked 
by a circle. Energies in units of /{kIb)- The number of 
translationally invariant states is much smaller than the total 
number of LLL states. 



to a description in terms of an excited rotating/vibrating 
quantal molecule. Specific examples of the reducibility of 
the JL and Moore-Read states to the molecular descrip- 
tion introduced in this paper are provided in Sections 
IIII CI and IIII Dl This is a surprising result, since these 
Jastrow based trial functions are widely described in the 
previous literature as being liquid-like in an essential way. 



B. Characteristic properties of the 
lowest-Landau-level spectra 



a given magnetic field B. This special form takes advan- 
tage of the simplifications at the limit of large B, i.e., 
when the relevant Hilbert space can be restricted to the 
lowest Landau level, given that fiujf) « ?iajc/2; the fre- 
quency ujo specifies the external harmonic confinement 
and LJc = eB / {mc) is the cyclotron frequency. Then the 
many-body hamiltonian reduces to 



e, global 
LLL 



N 



-a;2/4-^e/2)i + 5] 



K\Zi 



(1) 



where L = total angular momentum and 

z = X + iy. 

In the case of N rapidly rotating bosons (with ojq spec- 
ifying the external confinement of the two-dimensional 
harmonic trap and denoting the rotational frequency) , 



the corresponding Hamiltonian 
limit il/cjQ — > 1) is written as [46 1 




H 



b, global 

LLL — 



N 



NfiujQ + h{u}o - + ^ g5{z 



(in the 



(2) 



Since we will consider many-body energy eigenstates 
that are eigenstates of the total angular momentum as 
well, it follows that only the interaction terms are non- 
trivial in both Hamiltonians ([1]) and As a result, 
we will henceforth follow the practice of focussing on the 
simpler interaction-only LLL Hamiltonian 



N 



LLL 



(3) 



where v{zi — Zj) denotes the two-body interaction 
(Coulomb for electrons and repulsive contact potential 
for bosons) . The "ground states" of the Hamiltonian in 
Eq. dS]) coincide with the "yrast" band [36j . 

We remind the reader that the EXD method is based 
on the fact that the full LLL Hilbert space at a given 
total angular momentum L is spanned by the set of all 
possible uncorrelated permanents (for bosons) or Slater 
determinants (for electrons) made out from the Darwin- 
Fock zero-node single-particle levels (referred to also as 
orbitals) 



: exp(— zz*/2), 



(4) 



For completeness and clarity in the presentation, we 
briefly provide in this section a graphical illustration of 
some of the main characteristics of the LLL spectra, cal- 
culated via the exact-diagonalization approach. 

First we describe here the special form [ll|, [l^, |3^ 
of the many-body Hamiltonian used for calculating the 
global ground state of a finite number N of electrons at 



with li > 0. The position variable z is given in units of 
A = y'h/Jmujo) in the case of a rotating harmonic trap 
(with lateral confinement frequency ujo) or IbV^ in the 
case of an applied magnetic field B, with Ib — \/h/{mujc) 
being the magnetic length and uJc the cyclotron frequency 
[46j . (For details concerning the EXD method, see, e.g., 
Refs. [l^[ll|.) In the following, we use the convention 
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that an uncorrelated state is described by a single perma- 
nent (or Slater determinant) made out from the orbitals 
in Eq. ^ , which are characterized by good single-particle 
angular momenta li. 

A small part (sufBcicnt for our purposes here) of the 
EXD LLL spectra (as a function of L) are plotted in 
Fig. [T] for iV = 3 scalar bosons and in Fig. [2] for = 4 
spin-polarized electrons. As is usually done in the LLL, 
the one-body terms of the Hamiltonian (i.e., confining 
potential and kinetic-energy) were omitted [ill, [IE [13 , 
and the exact diagonalization involved only the two-body 
interaction [see Eq. 

In Figs. [T| and [2j the yrast bands ^ are denoted by 
a dashed line. Along the yrast bands there appear spe- 
cial cusp states denoted by a circle. The cusp states are 
important because they exhibit enhanced stability when 
the one-body terms of the Hamiltonian (i.e., external con- 
finement and kinetic energy) are added [see Eqs. (H]) and 
Qj. and thus they determine (47j the global ground states 
[lit ITsl f29l . lisj as a function of the applied magnetic 
field B or the rotational frequency of the t rap (for the 
correspondence between B and see Ref. [46[). In all 
studied cases [H [H [M [M [M Hll, Q (including both 
electrons and bosons up to A^ = 9 particles), the total an- 
gular momenta of the global ground states belong to the 
set of magic angular momenta given by Eq. (|30p below. 
We note that the emergence of these magic angular mo- 
menta are a direct signature of the molecular nature of 
the cusp states, a fact that further motivates our present 
investigations concerning the molecular description of the 
full LLL spectra beyond the electronic [l3, [ll, [3 cusp 
states. 

In the LLL, all many-body wave functions have the 
general formfH, M, 

W{z,,Z2,...,zn)\0), (5) 

where W{zi, Z2, ■ ■ ■ , zn) is an homogeneous polynomial 
of degree L (being antisymmetric for fermions and sym- 
metric for bosons). 

In Eq. ([5]), the symbol |0) stands for a product of Gaus- 
sians [see Eq. dU], i.e., 

N 

|0) =cxp(-^z.z*/2). (6) 

i=l 

To simplify the notation, this common factor will be 
omitted henceforth in the algebraic expressions and ma- 
nipulations [except in Eqs. (fT8| and ([33]) ]. where it is 
repeated for clarity]. Its contribution, however, is neces- 
sary when numerical results are calculated. 

A central property of the LLL spectra is the existence 
of a translationally invariant (TI) subspace [5, 25, 30] for 
a given L. This subspace is associated with a special 
subset of the general wave functions in Eq. ([S]), i.e., with 
wave functions having translationally invariant polyno- 
mials W{zi, Z2, ■ . ■ , z^). Specifically, the TI polynomials 
obey the relationship 

W{zi + C, Z2 + C, . . . , ZN + C) ^ W{zi, Z2, . . . , Zn), (7) 



for any arbitrary constant complex number c. 

The LLL states belonging to the TI subspaces are de- 
noted by gray solid dots (marked by an arrow; green 
online) in Figs. [T] and [5] The dimension D'^^{L) of the 
translational invariant subspace is much smaller than the 
dimension D^-^^{L) of the exact-diagonalization (EXD) 
[Tsl space (which is spanned by uncorrelated permanents 
or Slater determinants as discussed above). The remain- 
ing D^-^^{L) — D^^{L) states are spurious center-of-mass 
excitations, generated by multiplying the TI states with 
the operator z™, m = 0, 1, 2, . . ., where 

1 ^ 

i=l 

is the coordinate of the center of mass [49| . 

The energies of these spurious states coincide with 
those appearing at all the other smaller angular momenta 
5]. Thus (see TABLES H HI and|lV| 

i?TI(i) = Z?EXD(^)_^EXD(^_^)^ (9) 

We further note that for A^ particles (bosons or 
fermions) 

Dr{L)=Dy{L + N{N-l)/2), (10) 

where the subscripts b and / stand for bosons and 
fermions (electrons), respectively. N{N — l)/2 is the 
smallest value of the total angular momentum for spin 
polarized fermions in the LLL. 

C. Plan of the paper 

The paper is organized as follows: 

The analytic trial functions associated with pure ro- 
tations of bosonic molecules (i.e., the RBMs) are intro- 
duced in Section III Al followed by a description of the 
purely rotational electronic molecular functions (i.e., the 
REMs) in Section |IIBl 

Properties of the RBMs and REMs are discussed in 
Section HICl 

Section IIIDI introduces the general ro- vibrational trial 
functions (i.e., the RVMs). 

Case studies of the quantal molecular description of 
the LLL spectra are presented in Section IIIII In par- 
ticular. Section IIII Al discusses the case of A^ = 3 LLL 
scalar bosons, while Section IIII Bl discusses the case of 
A^ = 3 spin-polarized LLL electrons. The case of A^ = 4 
LLL electrons is presented in Section IIII CI along with 
an analysis of the Jastrow-Laughlin state (for fractional 
filling v — 1/3) from the viewpoint of the present molec- 
ular theory. Section IIII Dl investigates the case of A^ = 5 
LLL bosons, along with an analysis of the Moore-Read 
state according to the molecular picture. 

Section IIVI offers a summary and discussion. 

Finally, the Appendix discusses the rapid-convergence 
properties of the RVM basis. 
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We note that, going from iV = 3 to A'' = 5 particles, 
the molecular description requires consideration of suc- 
cessively larger numbers of isomeric molecular structures 
as elaborated in Section IIIII In particular, for = 3 
only one molecular isomer is needed, while three different 
molecular isomers are needed for N — 5. It is remark- 
able that these isomers are independent of the statistics 
(bosons or fermions). 



II. MOLECULAR TRIAL FUNCTIONS 

The molecular trial functions introduced in this paper 
are derived with the help of a first-principles methodol- 
ogy of hierarchical successive approximations which con- 
verge to the exact solution of the many-body Schrodinger 
equation [l^. Specifically, this methodology is based on 
the theory of symmetry breaking at the mean-field level 
and of subsequent symmetry restoration via projection 
techniques [l5|. In this Section, we present (and/or re- 
view where appropriate) this derivation in some detail. 



A. Purely rotational bosonic trial functions 
(RBMs) 

RBM analytical wave functions in the LLL for N 
bosons in two-dimensional rotating traps can be derived 
following earlier analogous derivations for the case of elec- 
trons [12||. Our approach consists of two steps: 

(I) At the first step, one constructs a permanent (Slater 
determinant for fermions) 'i'^ {zi, . . . , z^) out of dis- 
placed single-particle states u(zj, Zj), j — 1, . . . ,N that 
represent scalar bosons localized at the positions Zj, 
with (omitting the particle indices) z — x + iy = and 
Z + iY = Re''''. 



*^[z] = perm(7W^[z] 
with the matrix A^^[z] being 



M''[z] 



U(Zi,Zi) ... u^ZN,Zi) 



u{zi,Zm) ... u{zm,Zn) 



(11) 



(12) 



For the permanent of a matrix, we follow here the def- 
inition in Ref. [5l| , that is, the permanent is an analog 
of a determinant where all the signs in the expansion by 
minors are taken as positive. This definition provides 
an unnormalized expression for the permanent. If a nor- 
malized expression is needed, one has to multiply with a 
normalization constant 



where {pi,p2, • ■ • ,Pm} denote the occupations (multi- 
plicities) of the orbitals, assuming that there are M dis- 
tinct orbitals in a given permanent {M < N) |52i] . 



In the LLL, one can specifically consider the limit when 
the confining potential can be neglected compared to the 
effect induced by the gauge field. The localized u{z, Z) 
single-particle states (referred to also as orbitals) are then 
taken to be displaced zero-node Darwin-Fock states with 
appropriate Peierls phases due to the presence of a per- 
pendicular magnetic field [see Eq. (1) in Ref. [l^ . or due 
to the rotation of the trap with angular frequency £7. 
Then, assu ming a symmetric gauge, the orbitals can be 
represented [H, [13, H^l by displaced Gaussian analytic 
functions, centered at different positions Zj = Xj + Yj 
according to the equilibrium configuration of N classical 
point charges [53, Ull arranged at the vertices of nested 
regular polygons (each Gaussian representing a localized 
particle). Such displaced Gaussians are written as 

u{z , Zj) = (l/V7r) 

X exp[-\z - Zj\'^/2]exp[-i{xYj - yXj)], (14) 

where the phase factor is due to the gauge invariance. 
z = x + iy, and all lengths are in units of A in the case of 
a rotating trap or IbV^ in the case of an applied magnetic 
field; see Section llBl 

The localized orbital u{z, Z) can be expanded in a 
series over the complete set of zero-node single-particle 
wave functions in Eq. One gets [see Appendix A in 
Ref. ^56,]] 



u{z,Z) = J2Ci{Z)Mz), 



1=0 



with 



CiiZ) = {Z*yexp{-ZZ*/2)/VV. 



(15) 



(16) 



for Z ^0. Naturally, Co(0) = 1 and Ci>o(0) = 0. 

For an A^-particle system, the bosons are situated at 
the apexes of r concentric regular polygons. The ensu- 
ing multi-ring structure is denoted by (ni, n2, n^) with 
= N. The position of the j-th electron on the 



q-th ring is given by 



Z. 



Zgexp[i27r(l - j)/ng], 1 < j < n,. 



(17) 



The single permanent ^^[2] represents a static boson 
molecule. Using Eq. dTSl) . one finds the following expan- 
sion (within a proportionality constant): 



E 

ii=0,...,ijv=0 



Cl^{Zi)Ci.2{Z2) ■ ■ ■ C;„(Zjv) 



X P{h,h,-jNm, (18) 



(13) where P{liM,-,l 



N) 



perm[zii , 23 , • • • , 



'N J' 



the el- 
ements of the permanent are the functions z^ , with 
Zj^, . . . , zj^ being the diagonal elements. The Z^'s 
(with 1 < A: < iV) in Eq. ([HI) are the Zfs of Eq. (fT7)) . 
but relabeled. 
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In Eq. (ITSl) . the common factor |0) represents the prod- 
uct of Gaussians defined in Eq. To simphfy the no- 
tation, this common factor is usuahy omitted. 

(II) Second step: In the following, we will continue with 
the details of the complete derivation in the simplest case 
of a single (0, N) ring. Thus we consider the special case 



-3)/N 



l<j<N, 



(19) 



where R is the radius of the single ring. 

The Slater permanent breaks the rotational 

symmetry and thus it is not an eigenstate of the total 
angular momentum hL = ^J2jLi h- However, one can 
restore [l^, EBl the rotational symmetry by applying onto 
vE'^[z] the projection operator 



^^^^ 



27r 



d'ye 



(20) 



where HL are the eigenvalues of the total angular mo- 
mentum. 

When applied onto ^^[z], the projection operator Vl 
acts as a Kronecker delta: from the unrestricted sum in 
Eq. (fT8|) it picks up only those terms having a given total 
angular momentum L (henceforth we drop the constant 
prefactor h when referring to angular momenta) . As a re- 
sult the projected wave function = p^\|/^ is written 



as (within a proportionality constant) 

iiH Hn=L 



Pjh, . . .Jn) ,(4,^i^ + ...+4,^Ij^) 



(21) 



with ^2n{j -\)/N. 

We further observe that it is advantageous to rewrite 
Eq. (|21|) by restricting the summation to the ordered ar- 
rangements /i < ^2 < • ■ • < ^Af, in which case we get 



h+h^ Hn=L 



E 



P{Iu...,In) 



0<ll<l2<-<l 

perm[e*'^i'i , e*'^^'^ , . . . , e**"'"] 

Pl\p2l...pMl 



(22) 



The second permanent in Eq. ([22]) can be shown [57| to 
be equal (within a proportionality constant) to a sum of 
cosine terms times a phase factor 



jTr(N~l)L/N 



(23) 



which is independent of the individual Ij's. 

The final result for the (0, N) REM wave function is 
(within a proportionality constant): 



ll+l2 + ... + lN=L 

$g™(0,iV)[z]= E Ctihj2,...jN)veTm[z[\4,...,z'j^], (24) 

Q<ll<l2<...<lN 



where the coefficients are given by different expressions (i) For even iV, one has 

for even or odd number of bosons N. 



E cos { [{N - 1%, + {N- 3)/„, + ... + /, 



X.m V/.+i, -•••-(^-3)V_„-(A^-1)^.«]^}, (25) 



where the summation runs over the permutations of (ii) With the notation K = N — 1, the corresponding 
the set of N indices {1,2,..., N}. coefficients for odd N are: 
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N 



cos [[Kl^,+{K -2% 



f M 

n 

\k=l 



Pk'- 



21 — 2/ 



a{K} 



(26) 



where X]<T{if} runs over all permutations of iV — 1 indices 
selected out from the set {1,2,..., N} (of N indices). 

In both Eqs. ^ and the index M (with 1 < 

M < N denotes the number of different single-particle 
angular momenta 's (j = 1, 2, . . . , M) in the ordered list 
{li, I2, ■ ■ ■ , In} and the p^'s are the multiplicities of each 
one of the different Ij 's [occupations of the corresponding 
single-particle orbitals ipi^{z)]. For example, for iV = 4 
and {h = 2, = 2, 13 — 2,l4 — 5}, one has M = 2 and 
Pi = 3, p2 — 1; for {li — OJ2 — OJ3 = 0,^4 = 0}, one 
has M = 1 and pi = 4; for {h = 1, ^2 = 2, ^3 = 3, ^4 = 9}, 



one has M — A and pi ~ P2 — Ps = Pi ^ ^■ 

We further note that for both Eqs. ([IS]) and (EH) the 
total number of distinct cosine terms is N\/2 [5^, with 
the division by 2 following from the symmetry properties 
of cosine, i.e., from cos(— x) = cos(x). 

The RBM expresion for an (rii, ^2) two- ring configura- 
tion (with N ^ ni + n2) can be derived following similar 
steps as in the derivation of the expressions for the multi- 
ring REMs j^]. If Li+L2 = L, the final two-ring RBM 
expression is 



$^B^(ni,n2)[z] = 

ll-\-l2-^----\-lni—L'lUzi+l-^hzi+2-\-----\-iN—L2 

X! X! Cb{li,l2, lm)Cb{lm+i,lm+2, Zjv)pcrm[2:'\ , . . . , zj^], 

o</i<(2<...<i,ii o<;„i+i<;„2+2 <...<;« 

(27) 



where Cbih.h, ■ ■ ■ and Cb{lm+i,ln-,+2, ■ ■ ■ An) are 
calculated by applying the single-ring expressions of Eqs. 
and (j^. 

Generalizations of expression (1271) to structures with 
a larger number r of concentric rings involve for each 
g-th ring (1 < g < r): (I) Consideration of a separate 
factor Ch(Z„^_i+i, /„^_i+2, ■ ■ ■ , 'n,_i+n,); (11) A restric- 
tion on the summation of the associated Uq angular mo- 
menta, i.e., + ^«,_i+2 + • ■ ■ + ^n,„i+n, = Lq, with 



Refs. [13, following similar steps as those in Section 
III Al A determinant needs to be used, however, instead 
of a permanent, to conform with the antisymmetrization 
properties of the electronic (fermionic) many-body wave 
function. If Li + L2 = L, the final two-ring REM expres- 
sion is 



B. Purely rotational fermionic trial functions 
(REMs) 

The REM expresion for any (ni, 712, . . . , n^) multi-ring 
configuration (with N — X]q=i ^<?) derived earlier in 
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Cfih,l2, ■ ■ ■ , lni)C f{lni + l, lni+2, • ■ • , ^jv)det[z'^, , . . . , zj^], 

0<ii<i2<...<i,ii<i,ii +!<...<;» 

(28) 



where the fermionic coefScients C/(Zi, ^2, . . . , Z„j ) and 
C f{lni+i,lni+2, ■ ■ ■ ,In) are calculated by applying to 
each one of them the single-rmg [(0, N)] expression 




(29) 



It is straighforward to generalize the two-ring REM 
expression in Eq. psp to more complicated or simpler 
[i.e., (0, iV) and (l,iV — 1)] configurations by (I) consid- 
ering a separate factor C/(/„^_i+i, ?„^_,+2, ■ • . ,/„,_i+„J 
for each qth ring; (II) restricting the summation of the 
associated Uq angular momenta, i.e., /„^_j+i -1- ln^_i+2 + 
■■■ + = Lq, with X;g=i Lq ^ L. 

Apart from the permanent being replaced by a deter- 
minant, we note two other main differences between the 
REM and REM expressions. That is, for electrons (1) 
M = in all instances (single occupancy) and (2) a 
product of sine terms in the coefficients C/(. . .) replaces 
the sum of cosine terms in the coefficients Cb(. . .). 



C. Properties of RBMs and REMs 



The analytic expressions for ^^^^{ni,n2, ■ ■ ■ ,nr)[z] 
(the index RXM standing for either REM or REM) de- 
scribe pure molecular rotations associated with magic an- 
gular momenta 



£ = £o + ^ nqkq 

9=1 



(30) 



with kq, q = l,...,r being nonnegative integers; Cq — 
N{N — I)/2 for electrons and £o = for bosons. 

A central property of these trial functions is that iden- 
tically 



$™(ni,n2,...,n.)[z] =0 
for both bosons and electrons when 



£ ^ £o + ^ nqkq 

9=1 



(31) 



(32) 



This selection rule follows directly from the point group 
symmetries of the (rii, n2, . . . , r^r) multi-ring polygonal 



configurations. Indeed under condition ([5^ the Cb{- ■ ■) 
and C/(. . .) coefficients are identicall y ze ro, as can be eas- 
ily checked using MATHEMATICA In other words, 
purely rotational states are allowed only for certain an- 
gular momenta that do not conflict with the intrinsic 
molecular point-group symmetries. 

The yrast states corresponding to magic £'s [Eq. (|30l) ] 
are associated with the special cusp states described pre- 
viously in Figs. [Hand [2] and in Section ITbI Furthermore, 
the enhanced stability associated with the cusp states 
(see Section II Bl) is obviously due to the selection rule 
described by Eqs. and ((5^ . 

An important property of the REM and REM trial 
functions is their translational invariance (in the sense 
described in Section ITB|) . 



D. General ro-vibrational trial functions (RVMs) 

The RVM functions that account for the general exci- 
tations of the rotating molecules have the form (within a 
normalization constant): 



$r^(ni, 



(33) 



The index RXM stands for either REM, i.e., a ro- 
tating electron molecule, or REM, i.e., a rotating 
boson molecule. The purely rotational functions 
^^^^{ni,n2, ■ ■ ■ ,nr) have been described in detail in 
earlier sections. The product in Eq. p3p combines ro- 
tations with vibrational excitations, the latter being de- 
noted by Q™, with A being an angular momentum; the 
superscript denotes raising to a power m. Eoth $1^^'^ 
and are homogeneous polynomials of the complex 
particle coordinates zi, Z2, . . . , zjv, of order C and Am, 
respectively. The total angular momentum L = C + Xm. 
(5™ is always symmetric in these variables; is anti- 

symmetric (symmetric) for fermions (bosons). |0) is the 
product of Gaussians defined in Eq. ([5]); this product of 
Gaussians is usually omitted. 

The vibrational excitations Q\ are given by the same 
expression for both bosons and electrons, namely, by the 
symmetric polynomials: 



(34) 



where Zr is the coordinate of the center of mass defined 
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TABLE I: LLL spectra of three spinless bosons interacting via a repulsive contact interaction gS{zi — Zj). 2nd column: 
Dimensions of the EXD and the nonspurious TI (in parenthesis) spaces (the EXD space is spanned by uncorrelated permanents 
of Darwin- Fock orbitals). 4th to 6th columns: Matrix elements [in units of g/{TvA^), A = ^y h/ (mujo)] of the contact interaction 
between the correlated RVM states {k,m} [see Eq. ((36])]. The total angular momentum L = 3fc + 2m. Last three columns: 
Energy eigenvalues from the RVM diagonalization of the associated matrix of dimension D^^{L). There is no nonspurious state 
with L = 1. The full EXD spectrum at a given L is constructed by including, in addition to the listed TI eigenvalues [_D^^(L) in 
number], all the energies associated with angular momenta smaller than L. An integer in square brackets indicates the energy 
ordering in the full EXD spectrum (including both spurious and TI states). Seven decimal digits are displayed, but the energy 
eigenvalues from the RVM diagonalization agree with the corresponding EXD""^^ ones within machine precision. 



L 


^EXDpTI) 


{k, m} 


Matrix elements 






Energy eigenvalues (RVM diag. or EXD^') 





1(1) 


{0,0} 


1.5000000 






1.5000000[1] 




2 


2(1) 


{0,1} 


0.7500000 






0.7500000[1] 




3 


3(1) 


{1,0} 


0.3750000 






0.3750000[1] 




4 


4(1) 


{0,2} 


0.5625000 






0.5625000[2] 




5 


5(1) 


{1,1} 


0.4687500 






0.4687500[2] 




6 


7(2) 


{2,0} 


0.0468750 


0.1482318 












{0,3} 


0.1482318 


0.4687500 




0.0000000[1] 


0.5156250[4] 


7 


8(1) 


{1,2} 


0.4921875 






0.4921875[4] 




8 


10(2) 


{2,1} 


0.0937500 


0.1960922 












{0,4} 


0.1960922 


0.4101562 




0.0000000 


0.5039062[6] 


12 


19(3) 


{4,0} 


7.3242187x10"^ 


1.0863572x10"^ 


1.5742811x10"^ 










{2,3} 


1.0863572x10"^ 


0.1611328 


0.2335036 










{0,6} 


1.5742811x10"^ 


0.2335036 


0.3383789 


0.0000000 


0.0000000 0.5002441[13] 



in Eq. ([8|) and A > 1 is an integer positive number. Vi- 
brational excitations of a similar form, i.e., 



functions 



N 



(35) 



(and certain other variants), have been used earlier to 
approximate part of the LLL spectra. Such earlier en- 
deavors provided valuable insights, but overall they re- 
mained inconclusive; for electrons over the maximum 
density droplet [with magic £ — Co], see Refs. Q and 
0; for electrons over the = 1/3 (i' = ^o/^) Jastrow- 
Laughlin trial function [with magic £ = 3£o], see Ref. 
1C|; and for bosons in the range < L < N, see Refs. 

The advantage of Q\ |50j (compared to Q\) is that it is 
translationally invariant (TI) [l^ , a property shared 
with both $2^^ and In the following, we will 

discuss illustrative cases, which will show that the RVM 
functions of Eq. ([55]) provide a correlated basis (RVM ba- 
sis) that spans the TI subspace 0, [25l [30j of nonspurious 
states in the LLL spectra. 



III. MOLECULAR DESCRIPTION OF LLL 
SPECTRA 

A. Three spinless bosons 

Only the (0, 3) molecular configuration and the dipolar 
A = 2 vibrations are at play (as checked numerically), 
i.e., the full TI spectra at any L are spanned by the wave 



^RBM 
^3fc 



(0,3)Q™^{fc,m}, 



(36) 



with fc, m = 0, 1, 2, . . ., and L = 3k + 2m; these states are 
always orthogonal. 

Following Eq. (1^^ , a simplified analytic expression for 
the (0,3) RBM can be derived, i.e.. 



RBM 



(0,3) 



li+l2+h=C 

E 

o<ii<;2<i3 



CbiliMM) Perm[zl\ 



where the coefficients Ch{- ■ ■) are given by: 



h J3] 

(37) 




V l<i<i<3 



COS 



27r(Z, - Ij) 



(38) 



where 1 < M < 3 denotes the number of different single- 
particle angular momenta in the triad (hyhAs) and the 
Pf^'s are the multiplicities of each one of these different 
angular momenta. 

TABLE U provides the systematics of the molecular 
description for the beginning (0 < L < 12) of the LLL 
spectrum. There are several cases when the TI subspace 
has dimension one and the exact solution (I)'^'^^'^' coincides 
with a single {k,m} state. For L = the exact solution 
coincides with $0^"^ — 1 (Q° = 1); this is the only case 
when an LLL state has a Gross-Pitaevskii form, i.e., it is 
a single (normalized) permanent [see Eq. (j33|) ] given by 
|0) as defined in Eq. 
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For L = 2, we found cx Q2, with the index [i] 

indicating the energy ordering in the full EXD spectrum 
(including both spurious and TI states). Since [see Eq. 

m] 

1-1 (^^^ 

this result agrees with the findings of Refs. [2^, [S^l con- 
cerning ground states of bosons in the range < L < A^. 



For L = 3, one finds $oxact ^ ^ 



RBM 



Since [see Eq. 



m] 



RBM 



CX (21 - Zc){z2 - Zc){z3 - Zc), 



(40) 



this result agrees again with the findings of Refs. |25l.l59|. 

For L = 5, the single nonspurious state is an excited 
one, $1=2^'=' OC $P'^Q2- 

For L = 6 (i' = 1/2), the ground-state is found to be 



^cxact ^ 



_ 160 RBM 1q3 

(Zl - Z2)^(Z1 - 2:3)^(22 - 23)^, 



(41) 



N=4e L=18 



. [1] 



Q 
X 





I2> 




pure (0,4) REM 
(£=18) 



pure (1,3) REM 
(£=18) 



(0,4) REM (£=14) - 
{X=2, m=2) 



FIG. 

L = 



3: (Color online) CPDs for A'^ = 4 LLL electrons with 
18 {v = 1/3). Top row: The three lowest-in-energy 
EXD'^' states (see TABLE [TT|. Bottom row: The RVM trial 
functions associated with the largest expansion coefficients 
(underlined, see TABLE [lll| of these three EXD"^' states in 
the correlated RVM basis. See the text for details. The solid 
dot denotes the fixed point ro. Distances in nm. 



i.e., the bosonic Jastrow-Laughlin function for = 1/2 is 
equivalent to an RBM state that incorporates vibrational 
correlations. 

For L > N{N - 1) (i.e., v < 1/2), the EXD yrast 
energies equal zero, and with increasing L the degener- 
acy of the zero-energy states for a given L increases. It 
is important that this nontrivial behavior is reproduced 
faithfully by the present method (see TABLE |l| . 



B. Three electrons 

Although unrecognized, the solution of the problem of 
three spin-polarized electrons in the LLL using molecu- 
lar trial functions was presented by Laughlin in Ref. [60j . 
Indeed, the main result of Ref. ^6^] [see Eq. (18) therein] 
were the following wave functions (we display the poly- 
nomial part only) 



fc, m) (X 

{Za + iZbf" 



{Za - iZb) 



3k 



2i 



(4 + 4) 



2\m 



where the three-particle Jacobi coordinates are 

Zi+ Z2+ 23 



1/2 



Zl + Z2 



Z3 



Zb = -^(^1 - Z2). 



(42) 



(43) 



(44) 



(45) 



Expression (|^^ is precisely of the form ^'f^'^Q™, as can 
be checked after transforming back to Cartesian coordi- 
nateszi, 22, and 23. Thus the wave functions \k,m) of 
Ref. [6^ describe both pure molecular rotations, as well 
as vibrational excitations, and they cover the transla- 
tionally invariant LLL subspace. We note that the pairs 
of indices {k, m} are universal and independent of the 
statistics, i.e., the same for both bosons [Eq. p6p ] and 
electrons [Eq. ([l^ ]. as can be explicitly seen through a 
comparison of TABLE [J here and TABLE I in Ref. [g^. 

We further note that Laughlin did not present molecu- 
lar trial functions for electrons with A^ > 3, or for bosons 
for any A^. This is done in the present paper. 



C. Four electrons 

For A^ = 4 spin-polarized electrons, one needs to con- 
sider two distinct molecular configurations, i.e., (0, 4) and 
(1, 3). Vibrations with A > 2 must also be considered. In 
this case the RVM states are not always orthogonal, and 
the Gram-Schmidt orthogonalization is implemented. 

Of particular interest is the L = 18 case {v = 1/3) 
which is considered as the prototype of quantum- 
liquid states. However, in this case we found (see TABLE 
HI]) that the exact TI solutions are linear superpositions 
of the following seven RVM states [involving both the 
(0,4) and (1,3) configurations]: 



|1) 



<f,REM 
^18 



(0,4), 



|2)=$«f^(0,4)Qi 



|3) = $«^^(0,4)g4, 



|5) 
|7) = 



^RBM, 



1,3), 



|4)-<i>«™(0,4)Q6, 



$?r(l,3)Q^ 



(46) 
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TABLE II: LLL spectra of four spin-polarized electrons interacting via the Coulomb repulsion e'' l(ji\zi — Zj\) . Second column: 
Dimensions of the full EXD and the nonspurious TI (in parenthesis) spaces (the EXD space is spanned by uncorrelated deter- 
minants of Darwin-Fock orbitals). Last three columns: Energy eigenvalues [in units of e^/(K/s)] from the diagonalization of 
the Coulomb interaction in the TI subspace spanned by the trial functions ^^^f^,{Q,^)Q^ and $^+3't(l, 3)Q5;^ (RVM digonal- 
ization). Third to sixth columns: the molecular configurations (ni,n2) and the quantum numbers fc, A and m are indicated 
within brackets. There is no nonspurious state with L = 7. The full EXD spectrum at a given L is constructed by including, 
in addition to the listed TI energy eigenvalues [Z)^'(L) in number], all the energies associated with angular momenta smaller 
than L. An integer in square brackets indicates the energy ordering in the full EXD spectrum (including both spurious and 
TI states), with [1] denoting an yrast state. Eight decimal digits are displayed, but the energy eigenvalues from the RVM 
diagonalization agree with the corresponding EXD""^^ ones within machine precision. 



L 




[(ni, n2){fc, A, 


m}] 




Energy eigenvalues (RVM diag. or 


EXD"^^) 


6 


1(1) 


[(0,4){0,A,0}] 






2.22725097(1] 






8 


2(1) 


[(0,4){0,2,1}] 






2.09240211 [1] 






9 


3(1) 


[(1,3){1,A,0}] 






1.93480798[1] 






10 


5(2) 


[(0,4){1,A,0}] 


[(0,4){0,2,2}] 




1.78508849(1] 


1.97809256(3] 




11 


6(1) 


[(1,3){1,2,1}] 






1.86157215(2] 






12 


9(3) 


[(0,4){1,2,1}] 


[(0,4){0,2,3}] 


[(1,3){2,A,0}] 


1.68518201(1] 


1.76757420(2] 


1.88068652(5] 


13 


11(2) 


[(1,3){1,2,2}] 


[(0,4){1,3,1}] 




1.64156849(1] 


1.79962234(5] 




14 


15(4) 


[(0,4){2,A,0}] 


[(0,4){1,2,2}] 


[(0,4){0,2,4}] 


1.50065835(1] 


1.63572496(2] 


1.72910626(5] 






[(1,3){2,2,1}] 






1.79894008(8] 






15 


18(3) 


[(1,3){3,A,0}] 


[(1,3){2,3,1}] 


[(1,3){1,3,2}] 


1.52704695(2] 


1.62342533(3] 


1.74810279(8] 


18 


34(7) 


[(0,4){3,A,0}] 


[(0,4){2,2,2}1 


[(0,4){1,2,4}] 


1.30572905(1] 


1.41507954(2] 


1.43427543(4] 






[(0,4){0,2,6}] 


[(1,3){4,A,0}] 


[(1,3){2,2,3}] 


1.50366728(8] 


1.56527615(11] 


1.63564655(15] 






[(1,3){3,3,1}] 






1.68994048(20] 







The expansion coefficients of the three lowest-in-energy 
EXD states (labelled [1], [2], [4]; see TABLEHH in this 
RVM basis are listed in TABLE Hill One sees that for 
each case, one component (underlined) dominates this 
expansion; this applies for both the yrast state (No. [1]) 
and the two excitations (Nos. [2] and [4]). To further 
illustrate this, we display in Fig. [3] the conditional prob- 
ability (pair correlation) distributions (CPDs), 

P(z, zo) = ($1 •^(■^^ - - ^o)l*)' (47) 

for these three EXD"'"'^ states (top row) and for the RVM 
functions (bottom row) corresponding to the dominant 
expansion coefficients. The similarity of the CPDs in 
each column is noticeable and demonstrates that the sin- 
gle RVM functions capture the essence of many-body cor- 



TABLE III: iV = 4 LLL electrons with L = 18: Expansion 
coefficients in the RVM basis (labelled by the ]i)'s) for the 
three lowest-in-energy EXD'^^ states (labelled (1], (2], (4]; see 
TABLE nil)- The 4th column gives the RVM expansion coef- 
ficients of the corresponding Jastrow-Laughlin expression. 



RVM 


EXD" (1] 


EXD"^^ (2] 


EXD'^^ [4] 


JL 


|1> 


0.9294 


-0.3430 


0.0903 


0.8403 


|2) 


-0.1188 


-0.0693 


0.8930 


-0.1086 


|3) 


0.0067 


0.0382 


-0.2596 


0.0076 


|4) 


0.0137 


0.0191 


-0.0968 


0.0395 


|5) 


0.2540 


0.8486 


0.1519 


0.4029 


|6) 


0.0211 


0.0283 


0.3097 


0.0616 


17) 


-0.2387 


-0.3935 


0.0877 


-0.3380 



relations in the EXD states. Full quantitative agree- 
ment (within machine precision) in total energies can be 
reached by taking into consideration all seven RVM ba- 
sis states [see Eq. (|46p ]. Naturally, a smaller number of 
RVM states yields intermediate degrees of high-quality 
quantitative agreement. 

The celebrated Jastrow-Laughlin ansatz , 

$JL[^]= ]J (z.-z,f^+\ (48) 

\<i<j<N 

has been given exclusively an interpretation of a 
quantum- fluid state 0, [ll|. However, since the RVM 
functions span the TI subspace, it follows that any TI 
trial function (including the JL ansatz above and the 
compact CF states) can be expanded in the RVM basis. 
As an example, we give in TABLE IIIII (4th column) the 
RVM expansion of the JL state for iV = 4e and L = 18. 
One sees that, compared to the EXD yrast state (1st col- 
umn), the relative weight of the pure (0,4) REM (denoted 
by (1)) is reduced, while the weights of higher- in-energy 
vibrational excitations are enhanced. In this context, 
the liquid characteristics are due to the stronger weight 
of molecular vibrations which diminish the rigidity of the 
molecule. 

Of great interest also is the L = 30 [v = 1/5) case, 
which in the composite-fermion picture was found to be 
susceptible to a competition (17] between crystalline and 
liquid orders. However, we found that the exact nonspu- 
rious states for L = 30 are actually linear superpositions 
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TABLE IV: LLL spectra of five spinless bosons interacting via a repulsive contact interaction g5{zi — Zj) . Second column: 
Dimensions of the full EXD and the nonspurious TI (in parenthesis) spaces (the EXD space is spanned by uncorrelated 
permanents of Darwin- Fock orbitals). Last three columns: Energy eigenvalues [in units of g/{TvA'^)] from the diagonalization of 
the contact interaction in the TI subspace spanned by the trial functions $^^(0,5)^^, '^4k^hA)Q'x, and $1^^(2,3)0^ 
(RVM digonalization). Third to sixth columns: the molecular configurations (ni, 712) and the quantum numbers k, X and m are 
indicated within brackets. There is no nonspurious state with L = 1. The full EXD spectrum at a given L is constructed by 
including, in addition to the listed TI energy eigenvalues [D"^^ (L) in number] , all the energies associated with angular momenta 
smaller than L. An integer in square brackets indicates the energy ordering in the full EXD spectrum [including both spurious 
and TI states (EXD^^)], with [1] denoting an yrast state. Eight decimal digits are displayed, but the energy eigenvalues from 
the RVM diagonalization agree with the corresponding EXD^^ ones within machine precision. 



L 


D™(7?") 


[(ni, n2){k, X 


mil 




Energy eigenvalues (RVM diag. or 


EXD"^^) 





1(1) 


[(0,5)10, A,0ll 






O.UUUUUUUU[iJ 






2 


2(1) 


(0,5)10,2,111 






O. ( OUUUUUU[lJ 






3 


3(1) 


(0,5)10,3,111 






0. 1Z0UUUUU[1J 






4 


5(2) 


(1,4)11, A,0|1 


[(0,5)10,2,211 




9 '^nnnnnnn 11 


. 1 ( OUUUU [OJ 




5 


7(2) 


[(0,5){1,A,0}] 


[(0,5){0,5,1}] 




1.87500000[1] 


3.03125000(3] 




6 


10(3) 


[(1,4){1,2,1}1 


[(0,5){0,2,3}1 


[(0,5){0,3,2}] 


1.90664171(2] 


2.42925914(3] 


3.02347415(5] 


7 


13(3) 


[(1,4){1,3,1}1 


[(0,5){1,2,1}] 


[(2,3){1,2,2}] 


1.77354877[1] 


2.07258062(4] 


3.00543311(7] 


8 


18(5) 


[(1,4){2,A,0}] 


[(1,4){1,2,2}] 


[(0,5){1,3,1}] 


1.18821986[1] 


1.60795253(2] 


1.91688009(6] 






[(0,5){0,2,4}] 


[(2,3){2,2,1}1 




2.24905854[8] 


3.00273273(11] 




9 


23(5) 


[(1,4){1,5,1}1 


[(0,5){1,2,2}] 


[(2,3){3,A,0}] 


1.12895814(1] 


1.52195553(3] 


1.89102917(7] 






[(2,3){2,3,1}] 


[(2,3){1,3,2}1 




2.06074601(10] 


3.00082677(15] 




10 


30(7) 


[(0,5){2,A,0}] 


[(0,5){1,5,1}] 


[(0,5){0,5,2}] 


0.90026059(1] 


1.29362646(4] 


1.51398054(5] 






[(1,4){2,2,1}] 


[(1,4){1,2,3}] 


[(1,4){1,3,2}] 


1.66194766(8] 


1.95923264(14] 


2.12274862(17] 






[(2,3){2,2,2}] 






3.00035191(21] 






11 


37(7) 


[(0,5){1,3,2}] 


[(0,5){1,2,3}] 


[(1,4){2,3,1}] 


1.03755324(2] 


1.07552423(3] 


1.50429489(7] 






[(2,3){3,2,1}] 


[(2,3){2,5,1}] 


[(2,3){1,2,4}] 


1.58737738(10] 


1.94750687(18] 


2.03365831(20] 






[(2,3){1,4,2}] 






3.00012024(27] 
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47(10) 


[(0,5){2,2,1}] 


[(1,4){3,A,0}] 


[(1,4){2,2,2}] 


0.61480761(1] 


0.93028069(3] 


1.05066256(5] 






[(1,4){2,4,1}] 


[(1,4){1,4,2}1 


[(1,4){1,2,4}1 


1.34022509(10] 


1.50000444(11] 


1.50755634(13] 






[(2,3){4,A,0}] 


[(2,3){3,3,1}] 


[(2,3){1,3,3}] 


1.64279523(18] 


1.98164620(27] 


2.05477689(29] 






[(2,3){0,4,3}] 






3.00004768(36] 







of the following 19 [= D'^^L = 30)] RVM functions: 

$^™(0,4)Qi2-2/c^ ^ith fc = 0,l,2,3,4,5,6; 
*^-f3''.(l,3)Q^''^^ with fc = 2,4,6; 
$f^™(0,4)Q^^^/^ with fc = 0,3; 
<i>«f3^(l,3)Q^^ with fc = 2, 3, 4, 5, 6, 7, 8. (49) 



Diagonalization of the Coulomb interaction in the 
above TI subspace yielded an energy 0.25084902 /{kIb) 
per electron for the yrast state; this value agrees again, 
within machine precision, with the EXD result. The 
most sophisticated variants of the composite- fermion the- 
ory [including composite- fermion diagonalization (CFD), 
composite- fermion crystal (CFC), and mixed liquid-CFC 
states [ill, [13413] f^ll short in this respect. Indeed the 
following higher energies were found 17, 61j: 0.250863(6) 
(CFD), 0.25094(4) (mixed), 0.25101(4) (CFCL The CFD 
basis is not translationally invariant [1^, [l^- Conse- 
quently, to achieve machine-precision accuracy, the CFD 
will have to be performed in the larger space of dimension 
L»EXD(2, = 30) = 169. 



D. Five bosons 

For N = 5 spinless bosons, one needs to consider rovi- 
brational states [see Eq. p3p ] for three distinct molecular 
configurations, i.e., ^f^^{0,5)Q'^, $^^^(1, 4)g';^", and 
<^>fk^{2, 3)Q^. Vibrations with A > 2 must also be con- 
sidered; since the RVM states are not always orthogonal, 
a Gram-Schmidt orthogonalization is implemented. TA- 
BLE |IV] summarizes the quantal molecular description in 
the start of the LLL spectrum (0 < L < 12). 

Of particular interest is the L = 8 case, since it corre- 
sponds to the bosonic Moore-Read state, given [U, [3lj 
by the analytic expression: 

<i>^^i^[z]=5 n (^.-^if n (^fc-^o', (50) 

i<j£A k<leB 

where sets A and B contain [N - l)/2 and (iV + l)/2 
particles, respectively, if iV is odd. (For N even, both sets 
contain N/2 particles.) S symmetrizes over all possible 
ways of carrying such a division of N particles into the 
two sets. 

Based on sizeable overlaps with the EXD wave func- 
tions [13, HH, the bosonic Moore- Read states [3 are 
thought to represent the yrast states in the lowest Lan- 
dau level with the same angular momentum, i.e., L = 
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N=5 bosons L=8 




FIG. 4; (Color online) TV-body correlation functions for 
N = 5 LLL bosons with total angular momentum L = 8. (a- 
b): The yrast EXD'^' state (labeled [1]; see TABLE [rv| . (c): 
The Moore- Read state given by Eq. (|5Up . (d): Correspond- 
ing radial single-particle densities for the yrast EXD"""' state 
(solid line; red online) and the Moore-Read state (dashed 
line; blue online). Solid dots denote the four fixed points 
placed at ro exp[j7ri/2], j = 1,2,3,4 (with ro — 1.3A) in (a), 
(0,0) and ro exp[j7ri/2], j = 1,2,3 (with ro = 1.3A) in (b), 
and roexp[2ji'7rj/5], j = 1,2,3,4 (with ro = 1.02A) in (c). 
Crosses denote the expected position of the fifth localized bo- 
son according to the classical (1,4) polygonal-ring configura- 
tion in (a-b) and the (0,5) single-ring configuration in (c). For 
the (1,4) configuration, there are two non-equivalent ways of 
choosing the four fixed points. In (c), a node, associated with 
an octupole vibration (see text), develops at the expected 
position of the fifth particle. In the radial single-particle den- 
sities in (d), note the maximum (local minimum) at the origin 
of the solid line (dashed line) in agreement with the underly- 
ing (1,4) [(0,5)] configuration. 



1)2/2 for odd TV and L = N{N - 2)/2 for even N 

However, for = 5 and L = 8, we found (see TABLE 
ITV|) that the exact TI solutions are linear superpositions 
of the foUowing five RVM states [involving both the (0,5), 
(1,4), and (2,3) configurations]: 

ii) = $«^^(i,4), \2)=^r'^{iA)Qi 

|3) = $r''(0,5)Q3, \4)^<^>^^^{0,5)Ql 

|5^^^RBM(2,3)Q2. (51) 

The expansion coefficients of the yrast EXD"""^ state 
(labelled [1]; see TABLE |lVl in this RVM basis are listed 
in TABLE |Vl One sees that one component (labelled 
|1), underlined) dominates this expansion. To further 
illustrate this, we display in Fig.|3]the A^-body correlation 



functions for this EXD state, i.e., the quantity 

P{z; zl zl zl zl) (X |$(z; z?, z^, zg, zl)\\ (52) 

which gives the probability distribution of finding the 
fifth boson at a position z under the condition that the 
remaining four bosons are fixed at positions z^, j = 
1,2,3,4. The A^-body correlations exhibit an (1,4) con- 
figuration that corresponds to the dominant $^^'^(1,4) 
RVM component [component |1), defined in Eq. ([51]) ]. 
In particular, when the four fixed points are forming a 
complete square inscribed in a circle of radius ro = 1.3A 
[see Fig.lHa)], the fifth boson is localized around the ori- 
gin. When three fixed points coincide with only three 
vertices of this square, with the fourth being at the ori- 
gin, the fifth boson is localized around the fourth vertex 
of the square, as expected from the classical (1,4) config- 
uration. 

In contrast to the above, we found that the Moore- 
Read state [see Eq. ([50]) ] exhibits a drastically different 
behavior in its A^-body correlation function [portrayed in 
Fig. HJc)]. To analyze this behavior, we have expanded 
the Moore-Read state (for A^ = 5 bosons and L = 8) in 
the RVM basis of Eq. ([5T|) . As was the case with the 
Jastrow-Laughlin 1/3 function for A^ = 4 electrons, such 
an expansion is possible due to the fact that the MR state 
is translationally invariant. The corresponding expansion 
coefficients are listed in TABLE |Vl It is remarkable that 
one dominant coefficient (underlined) does appear, but 
it is associated with the RVM basis state $^^^(0,5)Q3 
(denoted as |3)), instead of the RVM state |1) that domi- 
nates the expansion of the EXD"^^ yrast state. This basis 
state 1 3) corresponds to a octupolar single-phonon vi- 
bration of a (0,5) polygonal configuration, and this is 
reflected in the A^-body correlation plotted in Fig. HJc). 
Indeed, with the four fixed points positioned at the ver- 
tices of a regular pentagon inscribed in a circle of radius 
ro = 1.02 A, the probability of finding the remaining bo- 
son is concentrated in the neighborhood of the fifth ver- 
tex that completes the pentagon (denoted by X), but in 
addition it exhibits a prominent node precisely at X . 

Fig.mjd) contrasts the radial single-particle densities of 
the EXD"""^ and MR states. Naturally, the radial single- 
particle densities provide a reduced amount of informa- 



TABLE V: A = 5 LLL bosons with L = 8: Expansion co- 
efficients in the RVM basis [Eq. ((5T|] (labelled by the |i)'s). 
The 2nd column gives the yrast EXD""^' state (labelled [1]; 
see TABLE [IVI). The 3rd column gives the RVM expansion 
coefficients of the corresponding Moore- Read expression [Eq. 

mi 



RVM 


EXD'^' [1] 


MOORE- READ 


|1> 


0.7879 


-0.5159 


|2> 


-0.1162 


0.1502 


|3> 


-0.6005 


0.7999 


|4> 


-0.0684 


0.1873 


|5> 


-0.0198 


-0.1908 
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tion regarding the correlation structure. However, note 
that they reflect the underlying (1,4) and (0,5) molecular 
configurations through the maximum at the origin of the 
solid line (EXD"^^ state) or the local minimum at the ori- 
gin of the dashed line (Moore- Read state), respectively. 

The above analysis provides a caveat against drawing 
conclusions by relying exclusively on overlaps (as is often 
the practice in the litterature of fast rotating ultracold 
bosons [llilHI)- Indeed, we can conclude that the Moore- 
Read state examined here disagrees in an essential way 
with the EXD many-body wave functions. 

We stress that the EXD yrast states with the same 
angular momenta as the MR states exhibit correlations 
that conform with the molecular structures associated 
with the magic angular momenta defined in Eq. (1301) . 
In particular, L = 8 (corresponding to the MR state 
for = 5) is a magic angular momentum associated 
with a (1,4) configuration (i.e., 8 mod 4 = 0). Another 
example is i — 12 (corresponding to the MR state for 

= 6). In this latter case, the EXD A^-body correlation 
function for the yrast state was studied in Ref. [29j , where 
it was found (see in paricular Figs. 5 and 6 therein) that 
it conformed to a (0,6) polygonal-ring configuration, in 
agreement with the fact that 12 mod 6 = 0. 



IV. SUMMARY AND DISCUSSION 

The many-body Hilbert space corresponding to the 
translationally invariant part of the LLL spectra of small 
systems (whether fermions or bosons, and for both low 
and high angular momenta) is spanned by the RVM 
trial functions introduced in Eq. ([33]). The yrast and 
excited states for both short- and long-range interac- 
tions can always be expressed as linear superpositions of 
these RVM functions. Thus the nature of strong correla- 
tions in the LLL refiects the emergence of intrinsic point- 
group symmetries associated with rotations and vibra- 
tions of molecules formed through particle localization. 
We stress the validity of the molecular theory for low 
angular momenta, where " quantum- liquid" physical pic- 
tures d, [nl, HH have been thought to apply exclusively. 
Our analysis suggests that liquid-type pictures, associ- 
ated with translationally invariant trial functions (e.g., 
the Jastrow-Laughlin, compact composite-fermion, and 
Moore- Read functions), are reducible to a description in 
terms of an excited rotating/ vibrating quantal molecule. 

In addition to the above conceptual advances, from 
a computational point of view, the introduction of the 
RVM correlated basis is promising concerning future 
computational developments. Indeed, it has the potential 
for enabling controlled and stepwise improvements of the 
variational wave function, in analogy with previous expe- 
riences from correlated bases in other fields of many-body 
physics (39I - I4H . Such developments may enhance com- 
putational capabilities for systems with a larger number 
of particles than it is currently possible. 

The main body of the present paper consisted of 



TABLE VI: Matrix elements of the contact interaction [in 
units of (//(ttA^)] between the RVM basis states for N = 5 
LLL bosons with L — 8. The notation for the RVM functions 
is the same as in Eq. (|5H) . Only four decimal points are 
shown. 





il> 


|2> 


|3> 


|4) 


|5) 


Il> 


1.3964 


1.5813x10"^ 


0.2567 


7.7995x10"^ 


0.1358 


|2) 




2.4933 


-0.1896 


-0.2631 


-0.3711 


|3) 






1.5554 


0.0000 


0.1930 


|4) 








2.4609 


0.2492 


|5> 










2.0588 



two parts. The general theoretical background of our 
methodology was presented in Section [TTl while Section 
mil was devoted to case studies. In particular: 

The analytic trial functions associated with pure rota- 
tions of bosonic molecules (i.e., the RBMs) were derived 
in Section lll Al followed by a description of the purely ro- 
tational electronic molecular functions (i.e., the REMs) 
in Section HTB] Properties of the RBMs and REMs and 
their association with magic angular momenta were dis- 
cussed in Section III CI The general ro- vibrational trial 
functions (i.e., the RVMs) were introduced in Section 

Concerning illustrative examples of the quantal molec- 
ular description of the LLL spectra. Section IIII Al dis- 
cussed the case of V = 3 LLL scalar bosons, while Sec- 
tion IIII Bl investigated the case of V = 3 spin-polarized 
LLL electrons. The case of V = 4 LLL electrons was 
elaborated in Section IlII CI along with an analysis of the 
Jastrow-Laughlin state (for fractional filling i' — 1/3) 
from the viewpoint of the present molecular theory. Fi- 
nally, Section llll DI studied the case of V = 5 LLL bosons, 
along with an analysis of the Moore- Read state according 
to the molecular picture. It was shown that the intrinsic 
correlation structure of the Moore-Read state disagrees 
strongly with that of the EXD wave function [in spite of 
having a rather good overlap with the EXD state, calcu- 
lated by us to be 0.913 (see also \24s])]. 

The Appendix discussed in detail the rapid- 
convergence properties of the RVM basis. 

In the case of bosons in harmonic traps, mean-field 
vortex solutions of the Gross-Pitaevskii (CP) equation 
have been proven most useful for the interpretation of 
experiments in a variety of situations of Bose-Einstein 
condensates with very large N; see, e.g., Refs. (62l - [63 |. 
Consequently, the results presented in this paper, and 
previously [lil . [28l . [29j , pertaining to crystalline- type cor- 
relations in finite systems with a small number N of LLL 
bosons (as well as the earlier prediction of the develop- 
ment of molecular-cr yst alline patterns in non-rotating fi- 
nite boson systems [65|) are unexpected. While these 
results, are of intrinsic interest for small systems, one 
may also inquire about the size-dependent evolution of 
the properties of the system with increasing N. This 
question, and in particular the possibility of a transition 
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TABLE VII: RVM-diagonalization total energies [in units of g/(jvK^)], relative errors, and overlaps at each intermediate step 
for the case of the yrast state for N = 5 LLL bosons with L — S (see Sec. IIIIDI) . The notation for the RVM functions is the 
same as in Eq. (|5ip . The composition of each step is shown in the second column. The last row displays the corresponding 
quantities for the Moore-Read function [Eq. H50|) ]. 



Step Composition Energy Relative error (%) Overlap 

~ \T) 1.3964006 (RVM) ITS 0.788 

2 |1)+|3> 1.2071388 (RVM) 1.6 0.990 

3 |1)+|3) + |2> 1.1947104 (RVM) 0.500 0.997 

4 |l>+[3) + |2)+|4) 1.1884821 (RVM) 0.02 0.9998 

5 |1) + |3)+|2) + |4)+|5) 1.1882199 (RVM) 0.00 1.000 



MOORE-READ 1.2658991 6.5 0.913 



TABLE VIII: RVM-diagonalization total energies [in units of e^/(K/_g)], relative errors, and overlaps at each intermediate step 
for the case of the yrast state for A'^ = 4 LLL electrons with i = 18 (see Sec. IIIIC|I . The notation for the RVM functions is 
the same as in Eq. (|46p . The composition of each step is shown in the second column. The last row displays the corresponding 
quantities for the Laughlin function [Eq. (|48|) ]. 



Step Composition Energy Relative error (%) Overlap 

~ |T) 1.3217670 (RVM) r23 0929 

2 |1)+15) 1.3174550 (RVM) 0.90 0.961 

3 \l)+\5}+\7) 1.3081859 (RVM) 0.19 0.992 

4 |l)-h|5)-h|7)-h|2) 1.3059258 (RVM) 0.015 0.99965 

5 |l)-h|5)-h|7)-h|2)-h|6) 1.3058052 (RVM) 0.006 0.99988 

6 |l)-hl5)-H|7)-H|2)-h|6) + |4) 1.3057413 (RVM) 0.001 0.99997 

7 \l)+\5} + \7)+\2) + \6) + \4) + \3) 1.3057290 (RVM) 0.000 1.00000 



JASTROW-LAUGHLIN 1.3105953 0.37 0.979 



of LLL bosons to mean-field GP behavior for large N, 
is of high importance conceptually and practically, and 
we expect that it will be the focus of future experimental 
and theoretical studies . 
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Appendix A: Convergence in the RVM basis 
1. The example of A = 5 bosons 

In this section, we analyze in detail the convergence 
properties of the RVM diagonalization for the yrast state 
of A'' = 5 bosons with L = 8 (a case associated with a 
Moore-Read function) that was discussed in Sec. IIII Dl 

In Table I VII we display the coupling matrix elements 
for the contact interaction between the RVM basis func- 
tions given in Eq. (j5ip . The convergence properties in 
the RVM basis cannot be immediately seen from an in- 
spection of these coupling matrix elements; while most of 
the off-diagonal elements in TABLE IVII are smaller than 
the differences between the associated diagonal elements, 
a couple of them are indeed larger. 



The fast-convergence properties of the RVM basis can 
be seen through a tabulation of the intermediate RVM- 
diagonalization total energies as the size of the RVM basis 
increases in successive steps. In TABLE IVIIi in addition 
to these intermediate RVM energies, we also display the 
corresponding relative error (relative to the EXD result) 
and the corresponding overlap with the EXD wave func- 
tion for the yrast state (denoted by the index [1] in Sec. 

inroi) . 

We stress that convergence (as a function of the num- 
ber of the RVM basis functions used in the calculation) 
is seen from TABLE IVIII to be achieved rapidly (i.e., al- 
ready with the use of only two basis functions one obtains 
a relative error of 1.6% for the energy eigenvalue, and a 
99% overlap with the exact eigenfunction) . In particular, 
we note that at the second step the RVM wavefunction 
is already more accurate compared to the Moore-Read 
function which exhibits a relative error of 6.5% for the 
energy and an overlap of 91.3% [see TABLE|Vll]. 



2. The example of N = 4 electrons 

In this section, we analyze in detail the convergence 
properties of the RVM diagonalization for the yrast state 
of A^ = 4 electrons with L = 18 (a case associated with 
a Jastrow-Laughlin function) that was discussed in Sec. 

Unci 

As previously, the convergence properties of the RVM 
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basis can be seen through a tabulation of the interme- 
diate RVM-diagonahzation total energies as the size of 
the RVM basis increases in successive steps. In TABLE 
IVIIIl in addition to these intermediate RVM energies, we 
also display the corresponding relative error (relative to 
the EXD result) and the corresponding overlap with the 
EXD wave function for the yrast state (denoted by the 
index [1] in Sec. ImCl) . 

Convergence (as a function of the number of the RVM 
basis functions used in the calculation) is seen from TA- 



BLE |VlTIl to be achieved rapidly (i.e., already with the 
use of only two basis functions one obtains a relative er- 
ror of 0.90% for the energy eigenvalue, and a 99% overlap 
with the exact eigenfunction) . In particular, we note that 
at the third step the RVM wavefunction is already more 
accurate compared to the Laughlin function: indeed the 
latter exhibits a relative error of 0.37% for the energy 
and an overlap of 97.9%, compared to 0.19% and 99.2%, 
respectively, in the case of the former [see TABLE IVIIIj . 
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